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Abstract 

The triply heavy baryons have a rather diverse mass range. While some of them pos- 
sess considerable production rates at existing facilities, others need to be produced 
at future high energy colliders. Here we study the direct fragmentation production 
of the fl CC c and Slbbb baryons as the prototypes of triply heavy baryons at the hadron 
colliders with different y/s. We present and compare the transverse momentum dis- 
tributions of the differential cross sections, p^ in distributions of total cross sections 
and the integrated total cross sections of these states at the RHIC, the Tevatron 
Run II and the CERN LHC. 
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1 Introduction 



Study of hadron production and decay has always been interesting. Histori- 
cally it has served to illuminate both the collider physics and the fundamental 
theories of interactions specially in the strong and electroweak sectors. Re- 
cently heavy hadrons have received great attention. Their structure is pre- 
dicted by the constituent quark model and, wherever light quarks are absent, 
they are nicely treated within the framework of the effective field theory and 
the perturbative QCD [1]. 

Singly and doubly heavy baryons, A's and S's , exhibit interesting properties 
and due to the involvement of the light degrees of freedom, they have been 
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studied in special models [2]. Triply heavy baryons are the heaviest composite 
states predicted by the constituent quark model. They are the last generation 
of baryons within the standard model. Essentially they are the fl ccc , fl cc b, flcbb 
and flbbb baryons. It has become clear that while the fl ccc and flbbb should be 
fragments of a c and a b quark respectively, the fl cc b and fl c bb may be produced 
each in a c or in a b quark fragmentation. The fragmentation of all triply heavy 
baryons has recently been studied and their production has been estimated at 
the CERN LHC [3]. However the wide range of their fragmentation probabil- 
ities (10~ 4 - 1CT 7 ), higher cross section of charm production at the Tevatron 
and RHIC and finally different acceptance cuts apart from different y/s for 
different colliders, suggest their production to be explored in other colliders as 
well. For the matter of simplicity, in this work we have chosen to investigate 
the production of the lightest (fi ccc ) and the heaviest (flbbb) as prototypes of 
the triply heavy baryons at the RHIC, the Tevatron Run II and the CERN 
LHC colliders. 



2 Fragmentation functions 



To evaluate the cross section of Q ccc and flbbb baryons at hadron colliders in a 
factorized scheme, we need their fragmentation functions. Indeed, it is possible 
to describe the fragmentation of these states by a single fragmentation function 
for Q — > flQQQ which has already been calculated [4]. In our model we have 
considered the emission of two gluons by the heavy quark Q, each producing 
a QQ pair. The three heavy quarks thus obtained form the Qqqq bound 
state leaving the heavy anti-quarks to form the final state jet. The bound 
state is characterized by the baryon decay constant Jb- The fragmentation 
process in leading order is described in Figure 1. Kinematically we have let 
the original heavy quark to keep its original transverse momentum and have 
ignored the respective motion of the constituents within the bound state. This 
perturbative picture is evaluated at the scale \i = fi a . The scale ji which is a 
scale at which such calculations are possible, is in the order of total mass of 
all final state particles, namely 5m c for c — > fl ccc and 5m b for b — > flbbb states. 
The result is the following fragmentation function 
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Fig. 1. The lowest order Feynman diagram illustrating the fragmentation process of 
a heavy quark Q into a triply heavy baryon Qqqq with identical constituent flavors. 

+109916/ - 49912/ + 20649/], (1) 



where a s is the strong interaction coupling constant evaluated at the pair 
of vertices of each gluon in the Figure 1, fs is the baryon decay constant 
which is defined in a similar manner to the meson decay constant, J'm and 
Cf is the color factor of the baryon state formed in the fragmentation of the 
heavy quark. Moreover, here we have defined £ = {q T )/m 2 with qr being the 
transverse momentum of the initial heavy quark and m is the heavy quark 
mass. The two functions f(z) and g(z) are defined as 



The function f(z) is the contribution of the energy denominator emerging 
from the phase space integration and the function g{z) is due to the quark 
and gluon propagators. 

The inputs for the fragmentation function (1) are the quark mass, baryon 
decay constant and the color factor. We have set m = m c = 1.25 GeV and 
m = nib = 4.25 GeV. For the decay constant and the color factor we have 
taken / B =0.25 GeV and C F = 7/6 for both cases of the f2 ccc and VL bbb states. 
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3 Inclusive production of a VLqqq baryon 



Inclusive production of a Qqqq baryon state in a hadronic collision is fulfilled 
in certain stages. Collision of hadrons provide the production of the required 
parton or partons which eventually fragment into a Qqqq state. The theoret- 
ical evaluation of this process is only possible at sufficiently high transverse 
momentum and taking advantage of the parton model factorization. All this 
is possible at a scale much higher than the calculable scale of the fragmen- 
tation functions. Therefore the fragmentation functions which are calculated 
at fragmentation scale, are evolved up to a scale at which the convolution of 
parton distribution functions, bare cross section of the initiating heavy quark, 
and the fragmentation function is possible. In other words, for the pp collision 
we may write 



da 
dp T 



[pp^n QQQ (p T ) + X] =YlJ dxidxjdzfi/pixi, n)fji P (xj, n) 

x Wj Q(Pt/z) + X, n)D Q „)]. (3) 



Here fi/ p (xi, /i) and fj/ p (xj,ii) are the parton distribution functions for the 
initial partons % and j carrying fractions Xi and Xj of the total momentum in 
the protons, a is the heavy quark production cross section and Dq^^ qqq (z, //) 
represents the fragmentation of the produced heavy quark into a triply heavy 
baryon. Note that here the scale \i in the parton distribution functions, sub- 
process cross sections and the fragmentation functions are set to be equal. 
For the parton distribution functions we have employed the parameterization 
due to Martin- Roberts-Stirling (MRS) [5], and have included the heavy quark 
production cross section up to the order of [6]. The s/s and acceptance cuts 
for the colliding facilities used in this work appear in Table 1. 

The production rates should not depend on the choice of the scale \x if both the 
production of high-energy partons and the fragmentation functions in all or- 
ders in a s in the perturbation expansion are included. However, only the results 
for next leading order parton production cross section and the leading order 
fragmentation function are available. Therefore the results obtained form (3) 
will depend on the scale \x. It is usual to choose the transverse mass of the heavy 
quark as the central choice of scale defined by /j, R = \Jpt 2 (parton) +mg 2 . If 
such a choice or sometimes multiple of it happens to be less than the frag- 
mentation scale, then the larger of (//, /i ) is chosen as appropriate scale. We 
have used the following form of the Altarelli-Parisi equation [7] to evolve our 
fragmentation functions 
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Vq^ d Q^QQq(z, V) = j —Pq^q(z/v, fi)D Q ^ QQQ (y, fi). (4) 

Here Pq^q(x = z/y, jj) is the Altarelli-Parisi splitting function. Note that only 
the term Pq~>q is included in (4). The reason is that the quark Q is assumed 
to be heavy enough to make other contributions irrelevant. The boundary 
condition on the evolution equation (4) is the initial fragmentation function 
Dq^qqq(z, h) evaluated at the fragmentation scale fi = fi . It should also be 
mentioned that the evolution of the fragmentation functions also sums up the 
logarithms of the order of /Ir/ uiq in the fragmentation functions. 

To check the sensitivity of our results, we have examined the behavior of the 
differential cross sections of Q ccc and flbbb at the RHIC, the Tevatron Run II 
and the LHC for different scales ranging from the fragmentation scale up to 
the scales below the Z° boson mass. It is seen that the sensitivity decreases 
with increasing scale. Such a study have led us to select appropriate scales for 
our further investigations. In this way we were motivated to choose the scales 
of AjiR and 6/ir for the Q ccc and Qm states respectively. 



4 Results and discussion 



In summary we have considered Q ccc (the lightest) and (the heaviest) 
triply heavy baryons to evaluate their production rates at the hadron collid- 
ers. To accomplish this, we have employed the next leading order results for 
parton production cross sections. The fragmentation functions used here are 
calculated in leading order perturbation theory. They provide reliable frag- 
mentation probabilities for the triply heavy baryons. To evaluate the cross 
sections, the well known patron model factorization at high transverse mo- 
mentum is employed. This procedure allows the complicated mechanism of 
hadron production to be treated in a factorized manner. This is possible in 
a scale which is higher than the scale at which the fragmentation functions 
are calculated. The Altarelli-Parisi evolution equation relates these scales. In 
the evolution of the fragmentation functions we have included only the Pq^q 
splitting function. Here, our evaluation of cross sections for the triply heavy 
baryons is devoted to different hadron colliders. Each collider with detection 
system has restrictions on the measurements of the transverse momentum and 
the rapidity of the particles. The so called acceptance cuts for the colliders 
considered here appear in Table 1. 

First we present the transverse momentum, pt, distributions of the differential 
cross sections at different hadron colliders for Q ccc and flbbb- They appear in 
Figures 2 and 3. Clearly, the distributions are sensitive to the choice of /i. 
Our choice of fj, — Ajj, R for Vt ccc and /i = 6/ir for Q bbb is at the region of 
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Table 1 

The center of momentum energy (\fs) and the acceptance cuts for the colliding 
facilities used in this work. The rapidity is defined as y = \ log{ [E — pi) / \E + pl)} . 





RHIC 


Tevatron Run II 


CERN LHC 


[GeV] 


200 


1960 


14000 


^ ut [GeV] 


2 


6 


10 


y < 


3 


1 


1 



Table 2 

The total integrated cross section in pb for f2 ccc and ilbbb baryons at different hadron 
colliders. The decimal places may not be significant within the uncertainties in- 
volved. 



Process 






RHIC 


Tevatron Run II 


CERN LHC 


c > Q ccc 


(atp = 




2758.3 


383.0 


308.0 


b — > Qbbb 


(atfi = 


6/^r) 


3.3 


6.0 


8.4 



the scale with minimum sensitivity. Although the cross section for a given p? 
differs up to three orders of magnitude from one collider to the other, they 
are still comparable. The difference is seen to grow with increasing p T . It is 
more interesting in the case of fl ccc where the distributions seem to converge at 
sufficiently low p?- Another important feature about these distributions is the 
rather high cross section of Q ccc which is more striking in the case of RHIC 
where low p^s are available. Figures 4 and 5 show the total cross sections 
for production of Q ccc and Qbbb with transverse momentum above a minimum 
value p™ in . Note that in Figures 2 and 3 only the range p T > p™ in = p™* were 
considered. It seems that the above mentioned general features, hold in the 
case of p™ m distributions, unless for that the difference in the differential cross 
sections for a given pt is not as much as that of p™ n distributions. Another 
general feature of our distributions is that the fall off decreases with increasing 
\fs and also with increasing p™ n or pt- 

We have also calculated the integrated total cross sections. They appear in 
Table 2. Note that the cross section for £l bb b increases with increasing y/s. But 
this is not the case for Vt ccc . Not only the order reverses in this case, but the 
cross section at RHIC is nearly one order of magnitude higher. 

What can we say about other triply heavy baryons? It seems that the baryons 
with at least two c or at least two b quarks, apart from the magnitude of 
their cross sections, behave in a similar fashion [3]. Therefore we expect that 
the Q cc b state produced in c or b quark fragmentation to have similar distri- 
butions as Q ccc . Likewise the Q c bb state emerging from a c or a b quark will 
behave like flbbb state. It is also interesting to note that our choice of Q ccc and 
Vtbbb with fragmentation probabilities of about 2 x 1CT 5 and 6 x 1CT 7 are not 
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Fig. 2. The pt distributions of the differential cross sections in nb for £l ccc production 
at the RHIC, the Tevatron Run II and the CERN LHC hadron colliders at the scale 
of \i = 4fi R . 



the states with maximum and minimum fragmentation probabilities. In other 
words, the lightest /heaviest of triply heavy baryons does not mean the state 
with maximum/minimum fragmentation probabilities. Indeed the fragmenta- 
tion probabilities for b — > Q cc t, and c — > Vt^ possess the maximum and the 
minimum fragmentation probabilities of 2 x 10~ 4 and 10~ 7 respectively. 

At the end we would like to discuss the reliability of our results. We have 
employed a fragmentation function obtained in a particular model. Without 
referring to its ingredients, it is worth mentioning that the predictions of this 
model is consistent with similar results in the case of doubly heavy baryons 
evaluated by Doncheski et al in [2] and also experimental results of [9]. This 
kind of comparison is presented in [3]. Moreover, similar fragmentation func- 
tions predict consistent fragmentation probabilities for J/ip and T states [8]. 
The remaining uncertainties are related to the calculation of the cross sec- 
tions. Here the the parton densities and the scales are relevant. These kinds 
of uncertainties are well described in the literature. Therefore, we believe the 
only factors which may alter our results are the consideration of higher order 
fragmentation functions and the heavy quark production cross sections. 
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Fig. 3. The same as Fig.2 but in the case of fi&y, production at the scale of 6[iR. 




Fig. 4. The p™ m distribution of the total cross section in pb for £l ccc production at 
the RHIC, the Tevatron Run II and the CERN LHC hadron colliders at the scale 
of n = A^R. 
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Fig. 5. The same as Fig. 4 but for the case of at the scale of jjl = 6^r. 
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